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ABSTRACT 


The Atmospheric Imaging Assembly instrument on board the Solar Dynamics Observatory has directly imaged the 
fast-propagating magnetosonic waves (FMWSs) successively propagating outward along coronal magnetic funnels. 
In this study we perform a numerical investigation of the excitation of FMWs in the interchange reconnection 
scenario, with footpoint shearing flow being used to energize the system and drive the reconnection. The modeling 
results show that as a result of magnetic reconnection, the plasma in the current sheet is heated up by Joule dissipation 
to ~10 MK and is ejected rapidly, developing the hot outflows. Meanwhile, the current sheet is torn into plasmoids, 
which are shot quickly both upward and downward. When the plasmoids reach the outflow regions, they impact and 
collide with the ambient magnetic field there, which consecutively launches FMWs. The FMWs propagate outward 
divergently away from the impact regions, with a phase speed of the Alfvén speed of ~1000 km s™!. In the k—w 
diagram of the Fourier wave power, the FMWs display a broad frequency distribution with a straight ridge that 
represents the dispersion relation. With the WKB approximation, at the distance of 15 Mm from the wave source 
region, we estimate the energy flux of FMWs to be E ~ 7.0 x 10° erg cm? s^! , which is ~50 times smaller than the 
energy flux related to the tube-channeled reconnection outflow. These simulation results indicate that energetically 
and dynamically the outflow is far more important than the waves. 
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1. INTRODUCTION 


In the past decade, observational analysis and theoretical 
modeling have shown that the solar atmosphere supports the 
oscillation or propagation of various modes of magnetohydro- 
dynamic (MHD) waves, such as slow-mode waves (Ofman et al. 
1997; De Moortel et al. 2000; Wang et al. 2002, 2009, 2013; 
Ofman & Wang 2002; Marsh et al. 2003; Ofman & Wang 2008), 
Alfvén waves (Tomczyk et al. 2007; De Pontieu et al. 2007; 
Cirtain et al. 2007; Okamoto et al. 2007; Jess et al. 2009; He 
et al. 20092, 2009b; McIntosh et al. 2011), and fast-mode waves 
(Schrijver et al. 1999; Williams et al. 2001, 2002; Nakariakov 
et al. 2003; McLaughlin et al. 2009). However, as a result of 
instrumental limitations, fast-propagating magnetosonic waves 
(FMWs) were scarcely observed in the past. To detect the fast 
magnetosonic waves, the field of view (FOV) of the instrument 
should be comparable to the wavelengths of the waves, which 
are the characteristic size of magnetic structures (e.g., length 
of coronal loops), and the cadences should be shorter than the 
typical periods of the waves, which range from a few seconds 
to several minutes (Nakariakov & Verwichte 2005). 

On board the Solar Dynamics Observatory (SDO), the new 
Atmospheric Imaging Assembly (AIA) has a 41' x 41’ full-Sun 
FOV at 1^5 resolution and high cadences up to 10s, which are 
all in favor of detecting fast propagating features. Using its EUV 
observations, Liu et al. (2010) firstly noted the FMWs in the low 
corona. As a first detailed analysis, Liu et al. (2011) found that 
the FMWs are driven quasi-periodically at the base of the flaring 
region and propagate outward along open funnels with a phase 
speed of up to ~2200kms~!. They also found that the main 


period of the FMWs coincides with quasi-periodic pulsations of 
the flare, suggesting a common origin of these phenomena. In 
greater detail, Shen & Liu (2012) investigated an FMW event 
occurring on 2011 May 30, and they further demonstrated the 
common physical origin of the FMWs and the accompanying 
pulsation flare, which is possibly driven by quasi-periodic 
magnetic reconnection generated by several physical situations, 
such as the repeated generation of plasmoids, as well as their 
coalescence. By reanalyzing this event with radio observations, 
Yuan et al. (2013) found that the FMW pattern is comprised of 
distinct wave trains with varying periods and wavelengths, and 
revealed again that the generation of FMWs should be tightly 
related to energy release in magnetic reconnections. Recently, 
Shen et al. (2013) presented observations of an FMW in multiple 
wavelengths and its interaction with other magnetic structures. 

So far, great modeling efforts have been devoted to un- 
derstanding the physics of FMWs in the solar magnetoatmo- 
spheres (Cargill et al. 1997; Wu et al. 2001; Bo et al. 2002; 
Nakariakov et al. 2003; Bogdan et al. 2003; Hasan et al. 
2005; Khomenko & Collados 2006; Hasan & van Ballegooijen 
2008; Schmidt & Ofman 2010; Fedun et al. 2011; Ofman et al. 
2011; Vigeesh et al. 2012; Pascoe et al. 2013). To investigate 
slow and fast magnetoacoustic-gravity monochromatic wave 
propagation, transformation, and interaction in the solar pho- 
tosphere and chromosphere, Bogdan et al. (2003) presented 
two-dimensional (2D), isothermal MHD simulations with the 
magnetic configuration consisting of a large unipolar magnetic 
flux concentration surrounded on each side by two smaller con- 
centrations of opposite-polarity field. To study the processes 
taking place around the magnetic canopy in a sunspot-like 
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structure, Khomenko & Collados (2006) conducted similar nu- 
merical simulations of magnetoacoustic wave propagation from 
the photosphere to the low chromosphere. With a global three- 
dimensional (3D) MHD code, Schmidt & Ofman (2010) simu- 
lated a coronal mass ejection of Gibson & Low flux rope type 
(Gibson & Low 1998) to excite a global FMW. According to 
the observation results presented by Liu et al. (2011), Ofman 
et al. (2011) performed a 3D numerical simulation to study 
the FMW generation, propagation, and trapping in the dipole 
magnetic field and gravitationally stratified density at coronal 
temperature. Pascoe et al. (2013) also performed 2D numeri- 
cal simulations of impulsively generated fast magnetoacoustic 
perturbations to explore the dispersion of fast magnetoacoustic 
waves in waveguides. These simulations either used the im- 
posed motions at the bottom boundary as an exciter to study the 
propagation of the FMWs, or used the large-scale perturbation 
to study the global nature of the FMWs. 

However, recent observations from the AIA indicate that the 
FMWSs may originate from magnetic reconnections associated 
with the flare (Liu et al. 2011, 2012; Shen & Liu 2012; Yuan 
et al. 2013). This work aims to provide a numerical investigation 
of this scenario. Utilizing the developed MHD model (Yang 
et al. 2013a, 2013b), we simulate the interchange flare scenario 
with the footpoint shearing flow to energize the system and 
drive the reconnection. The simulation results show that the 
plasmoid ejections from the X-point could provoke the FMWs 
propagating outward. By a comprehensive analysis of these, 
we find that the modeled FMWs share many features with the 
observed ones. 


2. NUMERICAL MHD MODEL 


The details of the 2.5 dimensional (2.5D) numerical MHD 
model used here have been described in Yang et al. (2013a, 
2013b). In brief, the model solves the conservation of mass, 
momentum, and energy, together with the induction equation. In 
the energy equation we account for ohmic heating with a switch- 
on resistivity of the same form as in Yang et al. (2013a). We 
also synthesize the coronal emission from the model as it would 
be observed with AIA/SDO using the temperature response 
function (Boerner et al. 2012), following the procedure of Peter 
& Bingert (2012). 

The simulation region spans 0Mm < x < 80Mm in the 
horizontal dimension and 0Mm < y < 75 Mm in the vertical 
dimension. Initially, the domain is covered by a nonuniform 
grid in both the x and y dimensions. The grid spacing is 
ôx = dy = 25km for OMm < y < 4 Mm, ôx = ôy = 50km 
for 4Mm < y < 6Mm, and ôx = dy = 100 km for y > 6 Mm. 
During the computation, the smallest grid cells are set in the 
region with the strongest current, with dx = dy = 12.5kmasa 
result of adaptive mesh refinement used in the numerical model. 
To solve the equations, we use a splitting-based finite-volume 
numerical scheme. The fluid part is solved by the second-order 
Godunov-type central scheme and the magnetic part by the 
constrained transport approach. The explicit second-order total 
variation diminishing Runge-Kutta time stepping is applied in 
time integration. Details about the numerics can be found in 
Feng et al. (2011) and Zhang et al. (2011). 

For the initial conditions, the temperature is set as a hyperbolic 
tangent function as follows: 


- 1 Y — Ycor 
T(y) = Tchr + 5 Foor = Tow) tanh | ——— J +1 , (1) 


Wir 
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where Toor and Torr are the respective temperatures in the corona 
and chromosphere. We set Toor = 1 x 10° K, T4, = 2 x 10^ K, 
Wr = 1.5Mm, and Yer = 3.0Mm. On the basis of the 
given function of temperature T(y), the distributions of the 
pressure p(y) and density p(y) are calculated from hydrostatic 
equilibrium and the ideal gas equation (Yang et al. 2013a, 
2013b). The initial velocity is 0. The initial magnetic field is 
constructed via the vector potential for a uniform background 
field and an infinite series of line dipoles as described by 
Edmondson et al. (2010) and used by Lynch & Edmondson 
(2013). For the boundary conditions, we use the same form 
as defined by Yang et al. (2013b), except for the footpoint 
shearing flow to energize the system. At the lower boundary 
of y = 0Mm, V; is set to 


0, I(x — xo) > d 
V,= 
lis sin(x (x — xo)/d)((x — xo —d?Y?, |(x — xo)| € d, 
Q) 


where xo = 27 Mm and d = 15 Mm. The maximum magnitude 
of the shear flows is (Vz)max = 15kms-!. We note that this 
value is close to the one employed by Antiochos et al. (1999) 
and smaller than the one used by Lynch & Edmondson (2013). 


3. NUMERICAL RESULTS 
3.1. Outflow from the Reconnection Site 


Figure 1 shows the resulting distributions of temperature T, 
number density N, horizontal velocity V. (see Animation 1), 
and vertical velocity V, (see Animation 2) for a snapshot at 
t = 7.13 minutes. From the figure, we can see that the shearing 
of the loop causes its field lines to bulge outward, compressing 
and deforming the null point. The shearing is in the z-direction, 
i.e., perpendicular to the plane shown in Figure 1. As a result, 
free magnetic energy and magnetic helicity are injected into 
the initial current-free system, and electric currents develop 
along the topological features, such as the null point and within 
the sheared loop. As the simulation progresses, the shearing 
motion elongates and strengthens the null-point current sheet, 
and magnetic reconnection with the creation of plasmoids due to 
tearing instability takes place. Then, the initially closed sheared- 
loop field lines open. Meanwhile, the initially open lines close 
down. Owing to Joule dissipation, the plasma near the null 
point is heated up to ~10 MK and is ejected to both sides of the 
null point, developing a pair of hot outflows, with one shooting 
upward and the other downward. 

From the distribution of the horizontal velocity V,, as well 
as Animation | and Animation 2 available in the online version 
of the journal, it can be seen that multiple arc-shaped wave 
trains are consecutively and alternately launched from the 
simulated flare. The waves propagate radically outward in nearly 
circular fronts and finally diminish in amplitude and leave the 
computation domain successively. It should be noted that this 
phenomenon has been shown, directly or indirectly, in some 
studies involving a vertical current sheet (Barta et al. 2008; 
Karpen et al. 2012), although those authors did not mention 
such waves explicitly. The distribution of the vertical velocity 
V, and the corresponding Animation 2 show that the slow waves 
also emanate from the reconnection site, just as in our simulation 
of a chromospheric jet (Yang et al. 2013b). 


3.2. Waves Launched from the Reconnection Site 


Figure 2 shows running difference images of temperature 
(top row), synthesized AIA emission at 171 A (middle row), and 
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Figure 1. Calculated distributions of temperature T, number density N, horizontal velocity Vy (see Animation 1), and vertical velocity V, (see Animation 2) 
at t = 7.13 minutes, where streamlines show the magnetic field lines. The white dashed box indicates the selected region for running difference image and 


Fourier analysis. 


(Animations of this figure are available.) 


synthesized AIA emission at 211 À (bottom row) in the white 
dashed box region outlined in Figure 1, which is at t — 7.63, 
7.83, and 8.03 minutes, respectively. To obtain the images at 
171 À and 211 A, we apply a point-spread function (Gaussian 
with 1/4 FWHM) to the synthesized images (076), rebind the 
images to the AIA pixel size, and finally plot the time series with 
the same cadence as AIA would observe. It can be seen that the 
arc-shaped wavefronts of V, are also observed in the images 
of temperature, the synthesized real emission at 171 À, and the 
synthesized real emission at 211 À. At t — 7.63 minutes, one 
arc-shaped train enters the selected region. At t = 7.83 minutes, 
it arrives at the center of view. At t = 8.03 minutes, the first 
arc-shaped train exits the view, and another one comes into view 
and behaves like the first one. 

Figure 3 shows temporal evolution of the perturbed veloc- 
ity components (x: solid; y: dotted), magnetic field compo- 
nents (same legend), and density at the point (x — 15Mm, 
y — 45Mm) on the path of the waves. The perturbed compo- 
nents are obtained after background detrending. In our case at 
this point, the fast wave propagates in the upper-left direction, 
while the background magnetic field aligns quite closely to the 
y-axis. Polarity relations of fast MHD waves predict that per- 
turbed quantities dB, and dV, are in negative correlation, dB, and 
dV, are also in negative correlation, and the perturbed magnetic 
field is in positive correlation with the perturbed density (Marsch 
1986). For fast MHD waves with plasma beta well below 1 like 
in the solar corona, the perturbed velocity component perpen- 
dicular to the background magnetic field (Bo) is much larger 
than the parallel component. In our case, at the position of our 
investigation, it is clear that dV, (almost parallel to Bo) has obvi- 
ously stronger amplitude than dV, (almost perpendicular to Bo). 


These features are also presented in Figure 4, and thus it is 
verified that these arc-shaped wave trains, emanated from the 
reconnection site, are FMWs. 

To illustrate how the simulated FMWs are driven, Figure 4 
presents the distributions of current density J., vertical velocity 
V,, and horizontal velocity V, (see Animation 3) at £ = 6.20, 
6.35, and 6.50 minutes, which are zoomed into the area around 
the reconnection region. In this figure, we can see that the 
collision between the ejected plasmoids and the field in the 
outflow region yields the FMW wave trains as shown in 
Figures | and 2. At f = 6.20 minutes, as a result of tearing 
instability, one plasmoid appears in the current sheet. Pushed by 
the reconnection outflow, this plasmoid moves upward along 
the current sheet with a speed of about 500kms~!. On its 
way, the plasmoid expands itself. Att = 6.35 minutes, it reaches 
the outflow region and quickly collides and reconnects with 
the magnetic fields existing there, rapidly releasing its plasma 
with high thermal pressure. Meanwhile, velocity perturbations 
in V, and V, appear in the outflow region. Another plasmoid 
is generated in the current sheet and is ejected downward. At 
t — 6.50 minutes, the velocity perturbations evolve into the 
outward-propagating arc-shaped wave trains, whose wavefronts 
could be fitted with an arc of a circle with its center at the 
site of the outflow region. At the same time, the downward- 
moving plasmoid collides and reconnects with the magnetic 
field existing at the outflow region, producing the arc-shaped 
wave trains, which propagate along and across the closed loops. 
When other plasmoids come into the outflow region, the process 
mentioned above will be repeated. Animation 3 illustrates that 
the wave pulses are generated once a plasmoid collides with 
the outflow region. When traveling outward, the wave pulses 
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evolve into distinguishable wavefronts and form multiple wave 
trains. It is also seen that the wave trains disappear after the 
coalescence into the outflow region. 

To see whether the eigenmode oscillations exist in the loops or 
not, we plot in Figure 5 the temporal profiles of Vproj (the velocity 
projected onto the slice) at different locations along the loops, 
with each plot showing the temporal profile at an equidistant 
position of 0.7 Mm along the loops. The fluctuations in Vproj 
shown in this figure lack evidence of eigenmode oscillations, 
ie. the loop does not oscillate as a whole at fundamental 
or harmonic frequencies, because no identical frequency and 
relevant phase difference can be identified from the oscillations 
at various positions. At this point it is unknown why eigenmode 
oscillations are absent in our simulation. This might be due to 
strong local damping or expansion of the 2D closed loop. While 
the eigenmode oscillation of the closed field lines is absent, 
flare-like magnetic reconnection is found to release the energy 
in phase with the wave generation. 

Figure 6 shows the calculated distribution of Alfvén speed 
V4 att — 7.13 minutes, with streamlines representing magnetic 
field lines. The figure exhibits a rather complicated distribution 
of Alfvén speed V4, where the Alfvén speeds vary in a wide 
range from tens to thousands of kilometers per second. Several 
strong V4 gradients across the magnetic field can be identified at 
the separatrices between different magnetic topologies. These 
gradients are formed self-consistently with the evolution of 
plasma and the magnetic field, rather than prescribed at the 
beginning. It can be seen that the waveguides confined by the 
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Figure 2. Running difference images of temperature (top row), synthesized real emission at 171 A (middle row), and synthesized real emission at 211 A (bottom row) 
in the white dashed box region shown in Figure 1, which is at t = 7.63, 7.83, and 8.03 minutes, respectively. 


boundary with the large V4 gradient do not possess a narrow 
spatial distribution, but extend over large angles, which is 
different from the results of the previous MHD simulations 
such as Ofman et al. (2011) and Pascoe et al. (2013). The 
existence of broad waveguides to confine waves consequently 
allows the simulated waves to propagate nearly isotropically 
with wavefronts in arcs with large central angles, rather than 
only appearing in narrow funnels with limited angular extents 
as observed by SDO/AIA (Liu et al. 2011; Shen & Liu 2012). 
This figure also indicates that the wave dispersion may be 
induced by the inhomogeneity of the Alfvén speed adjacent 
to the wave source region. As the wavefronts travel farther, they 
become distinguishable from each other (see Animation 3) 
while the propagation speed increases, and they finally form 
the multiple wave trains. 


3.3. Energy Transport of the Waves and Outflow 


To analyze the kinematics of the FMWs, we extract the 
values of dV. / V4o, dV,/ Vao, dT / To, and dN/No along the 
direction of wave propagation defined by the right boundary 
of the dashed box region shown in Figure 1, and we stack the 
obtained profiles in time sequence. This gives the time-distance 
diagrams as displayed in Figure 7, showing the steep, recurrent 
stripes resulting from the arc-shaped FMW wavefronts. The 
slopes of these stripes correspond to the propagating speeds of 
the FMWs. By adopting a linear fit with the points shown as the 
white cross in this figure, we measure the slope of each stripe 
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Figure 3. Temporal evolution of the perturbed velocity components (x: solid; y: dotted), magnetic field components (same legend), and density at the point (x — 15 Mm, 


y — 45Mm) on the path of the waves. 


and find that the first stripe has a slope of 1057 + 25 kms! , 
the second a slope of 960 + 21 km s^, and the third a slope of 
919 + 16kms~!. This gives an average propagating speed of 
the wave trains (Vpn) of about 1000 km s-!, which is the Alfvén 
speed in the model. Based on the parameters of these FMWs 
and the WKB approximation, we can estimate their energy flux 
by E = 0.5p(dv)? Vpn (with the nearly equal phase and group 
velocities seen below). Taking the number density of 105 cm? 
and dv of 30km s^! at the distance of 15 Mm from the source, 
marked as the black cross in Figure 1, we obtain a wave energy 
flux E ^ 7.0 x 10° erg cm? s~!. This is ~50 times smaller than 
the energy flux related to the tube-confined reconnection outflow 
at the nearby location marked as the white cross in Figure 1, 
which has a speed of ~350 kms !. 

Because the FMWs are generated impulsively, their transient 
nature might be interpreted as propagating wave packets in 
observations (Liu et al. 2011; Shen & Liu 2012). To investigate 
the recurrent period behavior of the FMWs, in Figure 8(a) we 
show a k—w diagram of a wave power. Like the analysis of 
FMW observations, an obvious linear steep ridge that describes 
the dispersion relation of the FMWSs is visible in the simulated 
k—q diagram. This steep ridge can be fitted with a straight line 
passing through the origin point, which gives average phase 
(Vpn = V/k) and group (vg, = dv/dk) velocities of the Alfvén 
speed of ~1000kms~!. Figures 8(b)-(d) are the snapshots 
of the Fourier-filtered Vx images with a narrow Gaussian 
function centered at the peak frequencies of 0.05 s^!, 0.15 s7}, 
and 0.25 s^!, respectively. These highlight the corresponding 
FMW wave trains at different frequencies. Figure 9 plots the 


k-averaged power spectrum to frequency, from the distribution 
displayed in Figure 8. We can see a continuous power spectrum, 
which basically decreases with increasing frequency without 
showing evident peaks. 


4. SUMMARY AND DISCUSSION 


In this study we performed a numerical investigation of the 
excitation of FMWs in the interchange reconnection in the solar 
corona. The footpoint-shearing flow is used here to energize the 
system and drive the reconnection. We find that the FMWs are 
initiated by plasmoid ejections from the X-point and that the 
generated FMWs are comparable to the observations. 

The modeling results show that the shearing of the loop causes 
its field lines to bulge outward, compressing and deforming the 
null point. As a result, the null-point current sheet forms and 
magnetic reconnection between sheared-loop field lines and 
open funnel lines takes place, with the creation of plasmoids 
due to tearing instability. The plasma near the null point is 
heated up by Joule dissipation to ~10 MK and is ejected to 
both sides of the null point. As a result of the collision between 
the impacting plasmoids and the field in the outflow region, 
the multiple arc-shaped fast-mode magnetosonic wave trains 
are consecutively and alternately launched from the outflow 
region, quickly propagating outward. Also, there is a one-to-one 
correlation between the energy release and the wave generation. 
In contrast, there are no evident eigenmode oscillations in the 
closed loop, which if they were present could be responsible for 
the wave excitation. 
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Figure 4. Distributions of current density Jz, vertical velocity Vy, and horizontal velocity V, at £ = 6.20, 6.35, and 6.50 minutes, respectively, zoomed into the area 
around the reconnection region. In this figure, streamlines show magnetic field lines, and black arrows in the middle row denote the velocity. 


(An animation of this figure is available.) 
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Figure 5. Temporal profiles of Vproj (the velocity projected onto the slice) at different locations along the loops. Panels (a)-(f) correspond to slices a-f shown in 
Figure 4 and Animation 3, respectively. In each panel the horizontal axis represents time, the vertical one Vproj; to distinguish the plots, they are stacked with a uniform 
increment. Each plot shows the temporal profile at an equidistant position of 0.7 Mm along the loops. The lower curves correspond to a position closer to the inside of 


the loop. 
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with a straight line passing through the origin, which gives the in Figure 5 may be negligible. In addition, prior observations and 
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Figure 7. Running difference time-distance diagrams of horizontal velocity d V; / V4o, vertical velocity dV, / Vao, temperature dT / To, and number density d N / No 
(for time differences of 1 s). The spatial direction is along the direction of wave propagation defined by the right boundary of the dashed box in Figure 1. V4o, To, and 
No are equal to 1000 km s7!, 10$ K, and 108 cm, respectively. The white crosses show the actual points that are used as the linear fits, and the yellow dashed lines 
indicate the wave propagation speeds. 
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Figure 8. Fourier analysis of FMWs in the dashed-box region shown in Figure 1. (a) Fourier power (k—w diagram) of a horizontal velocity field V, from t = 6.5 minutes 
to t = 10.0 minutes. (b)-(d) Vy sample images, which are Fourier-filtered with a narrow Gaussian function centered at the peak frequencies of 0.05 $71, 0.15s-!, and 


0.25 s71, respectively. 
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Figure 9. k-averaged power spectrum as a function of frequency, from the 
distribution displayed in Figure 8. 


simulations exhibit evidence of such reflection and bidirectional 
propagation (Liu et al. 2011; Ofman et al. 2011). The wave 
contribution from the other Y-type null to the perturbation in the 
selected region is also negligible in terms of relative fluctuation 
magnitude. 

In this work, we note that the simulated waves propa- 
gate isotropically from the source, not appearing in funnels 


with narrow angular extents as mostly observed by SDO/AIA 
(Liu et al. 2011; Shen & Liu 2012). This difference between the 
simulated and observed waves may affect their energy, propa- 
gation, and even frequencies and presents a possibility for im- 
proving our simulations in the future. In addition, we do not 
reproduce the observed power peaks of the k—« ridge. As sug- 
gested by Liu et al. (2011) and Shen & Liu (2012), the power 
peaks may be due to three-minute chromospheric oscillations 
or the leakage of the pressure-driven p-mode oscillations from 
the photosphere into the low corona. In the future, we may in- 
corporate these effects to investigate their possible modulation 
on reconnection and also the excitation of the FMWs. Another 
issue may arise as a result of the restriction of this work to a 
2.5D simulation. Although it is a simple control of the processes 
at work, the loss of one degree of freedom may affect the en- 
ergetics of the system and the excision of the FMWs, which 
has been mentioned in our previous studies (Yang et al. 2013b). 
Our future works will include an attempt to extend this work 
to a fully 3D simulation, which will more self-consistently en- 
force magnetic reconnection between the sheared loop and the 
ambient fields. 
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